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Abstract 

We compute the predictions of intrinsic charm for deep-inelastic scattering at high-Q^ 
and compare to HERA data. With the inclusion of constraints from low-energy data, 
enhancements beyond the predictions of standard structure functions are very modest, 
but peaking in the leptoquark mass variable is present near 200 GeV. Ultimately, the 
ability of HERA to probe the intrinsic charm hypothesis could be very substantial. 


HERA deep-inelastic scattering data with significant statistics at high and high x is now 
becoming available i0. The cross section exhibits an excess beyond the predictions based on 
standard QCD distribution functions. However, it is possible that unexpected contributions 
to the quark distribution functions might explain at least part of this excess p . In this Letter, 
we pursue the implications for the HERA data of an intrinsic charm (IC) component of the 
proton m. We hnd that intrinsic charm naturally predicts a peak in the leptoquark mass 
distribution in the vicinity of 200 GeV, but that existing constraints on intrinsic charm limit 
the size of the effect. As a function of leptoquark mass, the enhancement from intrinsic charm 
relative to standard model (SM) predictions is at most of order 15% in the neutral current 
(NC) case but can be as large as 75% in e~^p charged current (CC) scattering. Expectations 
for e~p NC and CC scattering are also given. 

The possibility of an intrinsic charm component in the proton bound state has substantial 
theoretical and phenomenological motivation. However, a conclusive observation or severe limit 
has remained elusive. In the intrinsic charm model, there is a \uudcc) component of the proton 
in which the c and c have become completely integrated members of the bound state, travelling 
coherently with the light quarks. In the minimal model, the \uudcc) wave function is taken to be 
proportional to the inverse of the light-cone energy difference rn^ — Yjijni/^i) ^ XcXc/{xc + xf) 
where the Xi are the light-cone momentum fractions and the light quark and proton masses 
are neglected compared to the charm quark mass. After squaring, the intrinsic charm Fock 
state probability distribution takes the form 

pp, 

dxudx',j^dxddxcdxc {xc + xf)^ ^ * ’ 


( 1 ) 




with n = 2. However, more extreme dependence on the light-cone energy difference is entirely 
possible and, indeed, very natnral, given that the strong interactions will tend to bind the 
intrinsic charm qnarks to the light qnarks so that all move with essentially the same velocity. 
To exemplify this sitnation, we will also consider n = 8 in Eq. (|l|). 

The inclusive charm quark distribution, c{x ), is obtained by integrating Eq. (|^ over all the 
Xi except X = Xc- For n = 2 and n = 8, one hnds 


c(x) = ^N^x^ 


L3 


(1 — a;)(l -I- lOx -f -f 2x(l + x)hix 


(n = 2), 


( 2 ) 


c{x) = 


+ 


35 + 1155x - 1575x2 - 11375x^ - 2450x^ + 490x^ - 96x® + 14x^ - x® 


X (l443 + 7161x + 5201x2 + {840 + 5880x + 5880x2} 


(n = 8) (3) 


respectively. We adopt a normalization such that the intrinsic charm Fock state component 
of the proton has 1% probability, implying = 36 (W = 2888028) for n = 2 (n = 8). The 
anti-charm distribution, c(x), is identical to c(x). Thus, at leading order, the net contribution 
from intrinsic c and c to the deep-inelastic electromagnetic structure function is F 2 ^ (x) = 
8xc(x)/9. Next-to-leading order QCD evolution corrections to the IC contribution to F 2 (x) 
are computed following the formalism outlined in Ref. p. (See also Ref. |p.) They are 
significant but evolve very slowly with Q^. 

In Fig. 1^, we illustrate results for the electromagnetic structure function, F 2 (x) at = 
25000 GeV2. First, we give the standard next-to-leading order prediction for F 2 ^ obtained by 
employing the recent MRS96(R2)|]^ distribution functions computed via perturbative QCD 
evolution, without inclusion of the IC component. Also shown is the n = 2 IC component 
alone and the SM plus n = 2 IC sum.Q A very modest enhancement in the vicinity of x ~ 0.2 
is observed. Results for the n = 8 model for the intrinsic charm component and corresponding 
total F 2 (x) are also shown. A much more substantial enhancement is observed, peaking in 
the vicinity of x 0.3. 

The low-energy direct EMC measurement p|] of F^lx) provides some evidence for intrinsic 
charm. In Fig. the EMC data for V = 168 GeV is compared to several predictions for 
^2- 0 The hrst is the perturbative “extrinsic” charm prediction of Ref. [^. In this approach, 
charm production in deep-inelastic scattering is computed to next-to-leading order using all 
relevant partonic-level subprocesses contributing to 'y*p —>■ cX. (In our comparison to this 
data, we employ the Aj^ = 0.239 GeV CTEQ3 0 distribution set and use scale pi = + 

in all distribution function and evolution evaluations |^.) The most important aspect of the 
subprocess approach is that charm-quark mass effects near threshold are correctly incorporated. 
This is not true of the perturbative charm distributions contained in the CTEQ3 distribution 

^In our QCD evolution computations for the intrinsic charm component, we employ = 0.344 GeV, 
as appropriate for MRS96(R2), and rric = 1.5 GeV. We use scale Q in evaluating distribution functions and 
QCD-evolution corrections to F 2 ^ (x). 

^At the values probed by EMC, the Z* contribution to F 2 is negligible compared to the 7 * contribution. 
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set itself; these are obtained by effectively assuming that the charm quarks are massless. 
That correct inclusion of mass effects is critical in the EMC energy range is shown by the 
curve in Fig. corresponding to the CTEQ3 charm distribution function. The other curves 
in Fig. 1^ illustrate the result of adding intrinsic charm components with n = 2 and n = 8 
to the extrinsic charm prediction.^ In agreement with Ref. |p, the n = 2 addition to the 
extrinsic charm distribution hts the data quite nicely. The n = 8 model provides equally good 
agreement with the data. Even larger values of n would still be consistent with the EMC data, 
but substantially larger probability for the IC wave function component would not. 


ISR data on pp 
intrinsic charm 


K + x 


provides additional evidence for and constraints upon 


11 , 12 


The production of Ac baryons can be viewed as arising from three 
sources. First, there is gg + qq fusion production of a cc pair, followed by c —Ac. Second, in 
inelastic pp collisions the c quark in the IC proton wave function component^ can fragment to 
a Ac. Finally, also in inelastic collisions, the c quark in the IC component of the proton wave 
function can coalesce with a u and d quark to produce a Ac; the xp distribution of the Ac 
produced in this manner is obtained by integrating over all the parton momentum fractions in 
Eq. ([^) with the constraint xp = Xu + Xd + Xc- Following the procedures of Refs. |^, it is 
found that the shape of dN/dxp at large xp > 0.5 is dominated by the third, i.e. coalescence, 
term. This is illustrated in Fig. ||, where dN/dxp is plotted for (i) fusion, (ii) fusion plus n = 2 
IC and (hi) fusion plus n = 8 IC. In each case, the prediction is normalized to a{xF > 0.5).Q 
We observe that the dN/dxp data requires the presence of an intrinsic charm contribution, 
and that the n = 2 shape is preferred over the n = 8 shape. It is for this reason that we have 
chosen not to consider still higher n values in Eq. (|l|) for the IC model, despite the fact that 
higher n values would yield larger effects at HERA. 

In order to compare to HERA data, we have chosen to compute the ratio of intrinsic 
charm plus the standard model perturbative prediction from MRS96(R2) to the MRS96(R2) 
prediction alone (a) as a function of after integrating over 0.1 < y < 0.9 and (b) as a 
function of M after integrating over y > 0.4, where M = y/xs is the eq leptoquark mass. 
Experimental results for these ratios have been presented by HI [|^ , and closely related results 
have been presented by Zeus p|. These ratios appear in Figs. || and respectively. Results 
for the n = 2 and n = 8 models are presented for e'^p and e~p, NC and CC reactions. 

The dependence in Fig. ^ shows that the n = 2 and n = 8 IC models predict e~^p 
NC enhancements that grow to roughly 6% and 12%, respectively, at high Q^. Much larger 
enhancements are indicated in the present HERA data. If these large enhancements persist, 
they would have to have a primary source other than intrinsic charm. Further, the relatively 
small intrinsic charm component could probably never be isolated from the dominant source 
without charm tagging. In the alternative scenario, where the enhancements decline with 


^ In the n = 2 IC case, the best two-component EC-I-IC fit is obtained for (1.27±0.06)F|"®-|-(0.92±0.53)F|^. 
The EC coefficient, 1.27, is not inconsistent with the expected K factor from higher order corrections, while 
the IC coefficient is consistent with a 1% IC probability. For the n = 8 curves, we have employed the same 
coefficients. 

^It is assumed that there is no extrinsic or perturbative c quark component to the proton wave function at 
the low momentum transfers involved. 

^Comparisons based on absolute normalizations are less reliable. 
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increased statistics, it is clear that the sensitivity to intrinsic charm could ultimately become 
very substantial. 

Data from other reactions potentially provide further opportunity for probing intrinsic 
charm. The additional possibilities include e~p NC scattering and e^p CC scattering. The 
intrinsic charm contribution to both NC and CC reactions is the same for e~^p and e~p scatter¬ 
ing. The relative IC enhancement in e~p NC scattering is then much smaller than in e~^p NC 
scattering due to the increased size of the SM cross section (deriving from the change in sign 
of the F 3 contribution). In contrast, relative enhancements in e~^p CC deep-inelastic scattering 
will be much larger, exceeding 60% and 150% at > 30000 GeV^ for the n = 2 and n = 8 IC 
models, respectively. These larger enhancements occur because the SM e'^p CC cross section 
is mainly proportional to the down quark distribution and is, therefore, suppressed compared 
to the SM e~^p NC cross section, which is largely proportional to the much larger (at large x) 
up quark distribution. Very small enhancements are predicted in e~p CC scattering. (That 
the above systematics should apply has been pointed out in Ref. [^.) 

The (SM-|-IC)/SM ratios as a function of the leptoquark mass, M, are particularly re¬ 
vealing. As seen in Fig. ^ the peak at moderate x values in the intrinsic charm distribution 
function (see Fig. |lD results in a substantial peak in the (SM-|-IC)/SM ratio for e+p NC scat¬ 
tering in the vicinity of M ~ 200 GeV, i. e. the same general location as the peak observed 
for data/SM by HI. Of course, the IC prediction for the height of the peak [(SM-|-IC)/SM 
= 1.06 and 1.12 for the n = 2 and n = 8 IC models, respectively] is much smaller than that 
observed; data/SM ~ 8 at M ~ 200 GeV in the HI data. Increased statistics will decisively 
determine whether or not the observed peak at M ~ 200 GeV in data/SM has anything to do 
with intrinsic charm. Observations of other reactions could prove very valuable. For example, 
the peaking in (SM-|-IG)/SM as a function of M predicted by the IG models will be larger 
in e^p GG scattering than in e^p NG scattering; in the former case, the (SM-|-IG)/SM ratio 
reaches a maximum of 1.3 (1.7) in the n = 2 (n = 8 ) IG models for M ~ 200 GeV. 

As a final point of comparison, we have also evaluated the effect of adding an extra com¬ 
ponent, 5u{x) = 0 . 02(1 — a;)°'^, to the valence u quark distribution function, as considered in 
Ref. P|. Results for such an addition (after evolution) are also shown in Figs. || and |^. We 
observe that, in the absence of other new physics contributions, data/SM as a function of M 
will allow one to easily discriminate between this possibility and intrinsic charm. 

In conclusion, we find that intrinsic charm is unable to explain enhancements in the e^p 
neutral current cross section as large as those seen at HERA, although a small peak in the 
vicinity of leptoquark mass M ~ 200 GeV is a natural prediction. In the absence of other 
new physics, HERA data will ultimately provide a sharp test of the intrinsic charm picture, 
especially once the charm component of F 2 is directly isolated. 
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Figure 1: Predictions for F 2 ^ (x) at = 25000 GeV^: solid — SM next-to-leading order 
prediction using MRS96(R2) 0 distribution functions only; dots — the n = 2 IC component 
alone; dot-dashed — SM plus n = 2 IC; dashes — n = 8 IC only; dash-long-dash — SM plus 
n = 8 IC. 
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Figure 2: EMC data for F^ix) at 17 = 168 GeV compared to: (i) dotdash — the extrinsic charm 
prediction of Ref. |^; (ii) dots — the CTEQ3 perturbative prediction; (iii) solid — EC+IC 
prediction for n = 2; (iv) dashes — EC+IC prediction for n = 8. 
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Figure 3: The xp distribution for pp —> Ac + X. Data from Ref. [|^ is compared to: (i) 
dotdash — gg + qq ^ cc fusion followed by c —> Ac; (ii) solid — fusion plus n = 2 intrinsic 
charm contributions; (iii) dashes — fusion plus n = 8 IC contributions. A 1% probability for 
the IC component of the proton wave function is used to £x the IC cross section. In all three 
cases, the overall normalization is hxed by a{xF > 0.5). 
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Figure 4: Predictions of various models for [da/dQ‘^]/[da^^/dQ"^] at HERA center-of-mass 
energy after integrating over 0.1 < y < 0.9. Results are shown for scattering and both 
neutral current (NC) and charged current (CC) scattering. Curve legend: solid — scat¬ 
tering and n = 2 IC; dotdash — e~^p, n = 8 IC; dashes — e~p, n = 2 IC; dots e~p, n = 8 IC; 
dash-dash-dot — e+p, 5u{x) = 0.02(1 — a;)°-^; dash-dot-dot — e~p, 5u{x) = 0.02(1 — a;)°-^. 
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Figure 5: Predictions for [da/dM\/[da^^/dM] at HERA center-of-mass energy after integrat¬ 
ing over y > 0.4. Notation as for Fig. the dash-dash-dot curve in the NC window is not 
shown since it is nearly identical to the dash-dot-dot curve. 
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